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Abstract: Drought is an inevitable condition with negative impacts in the agricultural and climatic sectors, 
especially in developing countries. This study attempts to examine the spatial and temporal characteristics of 
drought and its trends in the Koshi River Basin (KRB) in Nepal, using the standardized precipitation 
evapotranspiration index (SPEI) over the period from 1987 to 2017. The Mann-Kendall test was used to 
explore the trends of the SPEI values. The study illustrated the increasing annual and seasonal drought 
trends in the KRB over the study period. Spatially, the hill region of the KRB showed substantial increasing 
drought trends at the annual and seasonal scales, especially in summer and winter. The mountain region also 
showed a significant increasing drought trend in winter. The drought characteristic analysis indicated that 
the maximum duration, intensity, and severity of drought events were observed in the KRB after 2000. The 
Terai region presented the highest drought frequency and intensity, while the hill region presented the 
longest maximum drought duration. Moreover, the spatial extent of drought showed a significant increasing 
trend in the hill region at the monthly (drought station proportion of 7.6%/10a in August), seasonal 
(drought station proportion of 7.2%/10a in summer), and annual (drought station proportion of 
6.7% /10a) scales. The findings of this study can assist local governments, planners, and project 
implementers in understanding drought and developing appropriate mitigation strategies to cope with its 
impacts. 
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1 Introduction 


Drought is a recurrent climate phenomenon that occurs due to a reduction in normal precipitation 
amount and/or delays in the rainy season in a region for a specific time (Wilhite and Glantz, 1985; 
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di Lena et al., 2014). This phenomenon can last for months to years and presents a larger spatial 
extent than any other climatic hazard (Dabanli, 2018). Droughts may cause a significant shortage of 
the water supply and reduce agricultural production. Wilhite and Glantz (1985) categorized drought 
as a meteorological, agricultural, hydrological, and socioeconomic phenomenon based on the 
precipitation deficit over time. The causal effects of drought can be recognized at spatial and 
temporal scales by the intensity, duration, frequency, and severity characteristics (Wilhite, 2005; 
Touma et al., 2015; WMO and GWP, 2016). Sheffield and Wood (2007) studied the increasing 
drought duration from 1950 to 2000 at different latitudes at a global scale. Similarly, Trenberth et 
al. (2011) reported that the magnitude, intensity, frequency, and spatial extent of drought have 
doubled since the 1970s in different parts of the world, including South and East Asia. The Hindu 
Kush Himalayan region is the largest mountain range in the world. The freshwater towers of South 
Asia and some parts of Southeast Asia provide water to more than one billion people (Bolch et al., 
2012). This region is highly susceptible to natural hazards that have significant impacts on 
livelihood (Sharma et al., 2019). Nepal is one of the hazard-prone countries of the Hindu Kush 
Himalayan region along with its two large neighboring countries, China and India (Chen et al., 
2013; Liu et al., 2015; Das et al., 2016). The transboundary river basins in Nepal, such as Koshi, 
Gandaki, and Karnali, which originate from the high Himalayas to the plain land, sustain the 
livelihood of people in the mountain, hill, and low lands, which are facing extreme climate events 
(Sharma et al., 2019). The South Asian monsoon climate together with its strong interaction with 
extreme topographic variations within the short south-north distance predominantly affects the 
precipitation pattern in Nepal (Karki et al., 2009). The horizontal extension of hill and mountain 
ranges results in moist conditions on the south- and west-facing slopes and rain shadow effects on 
the northern sides of the slopes (Sigdel and Ikeda, 2012). The Intergovernmental Panel on Climate 
Change (IPCC) reported that the increasing trend in global temperature has largely impacted the 
precipitation pattern in different parts of the world (IPCC, 2018). A study by Manton et al. (2001) 
observed an increasing number of annual hot days and warm nights and a decreasing number of 
rainy days in the Asian region. Several studies have shown evidence of erratic precipitation patterns 
(Shrestha et al., 2000) and increasing maximum temperatures (0.06°C/a—0.12°C/a) in Nepal 
(Shrestha et al., 1999). Agrawala et al. (2003) revealed that Nepal will present an increment in 
mean temperature of 1.2°C and 3.0°C by 2050 and 2100, respectively. The increasing temperature 
may increase evapotranspiration and reduce soil moisture, which accelerates the high-risk climatic 
hazards in Nepal river basins. This risk is serious, especially in associated with extreme climatic 
events, such as drought (Easterling et al., 2000; Nijssen et al., 2001; Adnan et al., 2017). Since 
drought occurrence is primarily related to precipitation (decreasing) and temperature (increasing) 
factors, the aspects of drought intensity, duration, and frequency at spatial and temporal scales vary 
with the magnitude of these changing factors (Burke et al., 2006; Wanders and Wada, 2015; Jamro 
et al., 2019; Wang et al., 2019). Therefore, it is essential to study drought at spatial and temporal 
scales of Nepal. 

Drought indices are often used to study and monitor drought phenomena. Several drought 
indices have been developed for drought studies that can define and study drought at spatial and 
temporal scales (Mishra and Singh, 2010). Some of these indices are the Palmer drought severity 
index (PDSI) (Palmer, 1965), rainfall anomaly index (van Rooy, 1965), crop moisture index 
(Palmer, 1968), surface water supply index (Shafer and Dezman, 1982), and standardized 
precipitation index (SPI) (McKee et al., 1993). The development of the PDSI is considered the 
milestone in studying drought and one of the commonly used indices. However, its major 
limitations are that it is unable to capture short-term drought and does not possess multiscalar 
features in assessing drought (Alley, 1984; Yang et al., 2016; Wang et al., 2019). To overcome these 
limitations, researchers widely used the SPI for drought studies (Livada and Assimakopoulos, 
2007; Pai et al., 2011) because it is simple to calculate and can quantify the precipitation deficit for 
multiple time scales that reflect the drought impact on water resource availability (WMO, 2012). 
Because the SPI is based only on precipitation data, it is especially suited for studying 
meteorological drought (Portela et al., 2017; Wang et al., 2019). Many studies have indicated that 
temperature and evapotranspiration are also influencing factors for drought but are not considered 
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in the SPI; therefore, this method is not suitable for studying drought in the context of climate 
change (WMO, 2012; Wang et al., 2016; Jamro et al., 2019). Thus, the standardized precipitation 
evapotranspiration index (SPEI) was proposed by Vicente-Serrano et al. (2010); this index 
incorporates the merits of the SPI and PDSI, such as multitemporal scales, climate water balance, 
and evapotranspiration, and can reflect the influence of climate change on drought 
(Vicente-Serrano et al., 2010; Begueria et al., 2014). Therefore, this study considered the SPEI for 
drought analysis. 

Drought studies are relatively recent in Nepal. Ghimire et al. (2010) conducted qualitative 
assessments of climate-induced drought vulnerability at the household level of hill agriculture in 
Nepal. Kafle (2014) studied spatial and temporal drought assessments from 1982 to 2012 in 
western Nepal using the standardized reconnaissance drought index. Adhikari (2018) reviewed 
climate-induced drought hazards in mid-hill farming systems in western Nepal. Baniya et al. (2019) 
studied spatial and temporal analyses of drought in Nepal for the period 1982-2015 using the SPI 
and vegetation condition index. Hamal et al. (2020) studied the impact of drought on crop yields 
across Nepal using the SPEI method. 

There are three major river basins in Nepal, namely, Gandaki, Karnali, and Koshi. Dahal et al. 
(2016) conducted a drought risk assessment in the Gandaki River Basin of central Nepal from 1981 
to 2012 using the SPI. Similarly, Khatiwada and Pandey (2019) characterized the 
hydrometeorological drought in the Karnali River Basin of western Nepal from 1981 to 2014. In 
this study, several drought indices were used and the results revealed that the SPI is correlated with 
other indices and is capable of capturing hydrometeorological drought, which indicates that in the 
Karnali River Basin, the significance of other drought factors is negligible; under such conditions, 
drought is mostly contributed by precipitation variations (Vicente-Serrano et al., 2012). 
Nevertheless, several studies have shown the role of warming-induced drought in the context of 
climate change (Li et al., 2020; Pei et al., 2020) and the impacts of drought on agriculture (Chen et 
al., 2016; Wang et al., 2016; Hamal et al., 2020). 

A wide range of water-related hazards have affected the Koshi River Basin (KRB), and among 
these hazards, drought has affected large areas (Chen et al., 2013). The river basin receives 
approximately 900.0—1200.0 mm of annual precipitation, and the annual mean evapotranspiration 
is approximately 1500.0—2000.0 mm (Chen et al., 2013). Therefore, this river basin always has a 
water shortage. The basin shows increasing weather extremes (Shrestha et al., 2017), and seasonal 
precipitation is unevenly distributed. Approximately 75% of the annual precipitation occurs during 
the summer monsoon (from June to September), 10% occurs during October and November, and 
15% occurs during the winter and spring seasons (Agarwal et al., 2014). The daily minimum and 
maximum temperatures are increased by 0.3°C/10a and 0.1°C/10a, respectively, and consecutive 
dry days are increasing (Shrestha et al., 2017). The monthly average maximum and minimum 
temperatures range from 27.0°C to 32.0°C and from 0.0°C to 5.0°C, respectively. The studies cited 
above also projected that the annual and seasonal minimum and maximum temperatures will have 
an increasing trend in the future in the entire basin. Approximately 70% of the population in the 
river basin relies on rainfed agriculture for their livelihood; therefore, the uneven precipitation 
distribution and increasing temperatures cause water stress in the basin (Dixit et al., 2009). 

Several studies have been conducted in the KRB regarding these variations in precipitation 
(Agarwal et al., 2014), temperature (Agarwal et al., 2016), water resources (Bharati et al., 2019), 
climate change, and agricultural production (Neupane et al., 2013; Bhatt et al., 2014). Bharati et al. 
(2014) mentioned the future impact of climate change and the inclusion of extreme climate events, 
such as drought, in the planning process. A study of Wu et al. (2019) was conducted in the entire 
basin (including the transboundary region) on drought using the satellite-based crop water shortage 
index from 2000 to 2014. Shrestha et al. (2017) estimated the SPI values in the KRB (Nepal's part) 
to compare the precipitation-based drought severity and assess the accuracy of two climate hazard 
group satellite precipitation estimates. As mentioned earlier, increasing temperature is also one of 
the factors affecting drought occurrence, and the KRB is already facing increasing weather 
extremes. In this scenario, studying drought in the context of climate change based on 
meteorological indices is also essential. Thus, this study aims to (1) examine the spatial and 
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temporal variations in drought in different ecological regions of the KRB from 1987 to 2017; and 
(2) analyze drought characteristics using the SPEI at a three-month time scale from 1987 to 2017. 
This study evaluated drought characteristics using the SPEI based on ground station meteorological 
data over 31 years in the KRB, Nepal. The investigation of the spatial and temporal variability of 
drought in this study can provide valuable information for decision-makers to devise plausible 
policies, make climate-smart decisions, and implement these policies and decisions. 


2 Materials and methods 


2.1 Study area 


The study area is located at 26°49’-27°59'N latitude and 86°55'-87°10’E longitude within the KRB 
(Koshi River Basin) in the central and eastern regions of Nepal. The KRB is a transboundary river 
system with a total area of 8.73x10* km? and covers five counties in China, 27 districts in Nepal, 
and 16 districts in India (Bharati et al., 2019). In Nepal, the basin covers three distinct ecological 
regions in general, namely, the mountains, hills, and Terai, from north to south, and the elevations 
range from 60 to 8849 m a.s.l. (Fig. 1) (Khanal et al., 2018). The climate within the KRB (Nepal) 
varies from tropical in the Terai to alpine on high mountains. The climate in the hills is subtropical 
warm temperate (Dixit et al., 2009). Generally, the temperature decreases from south to north, 
although the topography causes spatial variations. The annual mean temperature in the mountains, 
hills, and Terai is 16.0°C, 20.0°C, and 25.0°C, respectively. The annual precipitation decreases 
from 1931.0 to 536.0 mm. Among the 27 districts in the KRB (Nepal), we selected 16 districts for 
our study based on data availability. Among these districts, three are located in the mountains, 
seven are located in the hills, and six are located in the Terai ecological zone (Fig. 1). The major 
climatic features of the mountain, hill and Terai regions estimated from the available data (1987— 
2017) of these studied districts are illustrated in Table 1 (Norbu, 2004; Dixit et al., 2009). 
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Fig. 1 Overview of the Koshi River Basin (KRB), Nepal (a) as well as the distributions of weather stations in the 
KRB (b) 
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Table 1 Major climate features of the mountain, hill and Terai regions (derived from the selected districts) of the 
Koshi River Basin (KRB), Nepal 


: Annual mean minimum Annual mean maximum Annual precipitation ; Elevation 
Region Climate 
temperature (°C) temperature (°C) (mm) (m) 
Mountain 11.3 22.3 1953.5 Gool temperate to >2000 
alpine 
Hill 13.0 24.2 1641.4 Subtropical to wam = 30ġ-oggg 
temperate 
Terai 19.2 30.6 1637.2 Tropical <300 


2.2 Data sources 


The data employed in this study included the monthly precipitation and minimum and maximum 
temperatures. The data were collected from the Department of Hydrology and Meteorology, 
Government of Nepal. For the current study, we considered the longest period of climate data 
available from the maximum number of weather stations. Altogether, 31 years of data from 23 
weather stations between 1987 and 2017 were collected. The 23 weather stations are located within 
16 districts of the KRB (Nepal). Among these 23 weather stations, 3 are in the mountain region 
districts (hereafter mountain region), 11 are in the hill region districts (hereafter hill region), and 9 
are in the Terai region districts (hereafter Terai region). The details of weather stations are given in 
Table 2. For the consideration of basin- and regional-level distributions of drought, we calculated 
the average value of climate data at multiple time scales from all weather stations in the basin and 
respective regions (mountain, hill, and Terai) and obtained the SPEI values (Sigdel and Ikeda, 
2010; Liu et al., 2018; Hamal et al., 2020). We assumed that the collection of the longest period of 
data from the maximum available stations for the drought study can help to minimize the 
uncertainties regarding the spatial variation in drought in the region. 


Table 2 Details of the weather stations selected from different regions in the KRB, Nepal 


Region District Station name Elevation AMI CV-AMP (%) PEI CV-PET (%) 
m) (mm) (mm) 
Mountain Dolakha Jiri 1877 2441.4 10.3 1222.3 2.3 
Sankhuwasabha Chainpur 1277 1444.2 14.3 1367.8 4.7 
Taplejung Taplejung 1744 1974.9 13.2 1161.3 4.3 
Hill Kathmandu Panipokhari 1329 1504.9 13.2 1449.1 6.1 
Kathmandu Kathinaridu 1337 1480.4 14.8 1476.0 2.6 
Airport 
Lalitpur Khumaltar 1334 1163.4 16.7 1394.6 2.4 
Lalitpur Godavari 1527 1708.6 18.0 1273.1 6.3 
Bhaktapur Nagarkot 2147 1863.3 14.4 1129.7 9.3 
Makwanpur Hetauda 452 2365.3 17.9 1607.6 3.8 
Kavre Dhulikhel 1543 1483.1 18.8 1230.6 8.0 
Kavre Panchkhal 857 1121.9 17.6 1621.1 2.8 
Okhaldhunga Okhaldhunga 1731 1781.4 11.4 1185.4 7.3 
Dhankuta Pakhribas 1720 1530.5 12.4 1110.1 2.5 
Dhankuta Dhankuta 1192 967.0 19.9 1261.1 9.8 
Terai Bara Simara 137 1814.7 22.0 1609.1 2.9 
Bara Parwanipur 87 1591.9 21.8 1599.3 3.5 
Sarlahi Karmaiya 139 1834.0 20.3 1549.3 2.8 
Sarlahi Manusmara 90 1438.8 30.2 1608.5 2.8 
Dhanusa Janakpur 76 1523.2 26.8 1540.5 2.7 
Siraha Lahan 110 1368.1 22.5 1485.2 12.4 
Saptari Rajbiraj 68 1500.0 26.4 1601.8 55 
Saptari Phattepur 83 1731.6 23.9 1693.5 4.5 
Sunsari Tarahara 120 1976.1 17.6 1553.3 2.7 


Note: AMP, annual mean precipitation; CV, coefficient of variation; PET, potential evapotranspiration. 
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2.3 Methods 


In this study, we used the SPEI index and Mann-Kendall (M-K) test to quantify drought and detect 
its trend in the KRB. The SPEI is considered an improved drought index to study drought 
conditions (Begueria et al., 2014). Similarly, the M-K test is a widely used nonparametric approach 
to detect trends in time series data (Salmi et al., 2002). The details of these methods are illustrated 
in the following sections. 


2.3.1 Standardized potential evapotranspiration index (SPEI) 


The SPEI is a drought index developed by Vicente-Serrano et al. (2010) that uses the precipitation 
and temperature components to calculate the index. The strength of this index is the inclusion of 
temperature components along with precipitation to account for its effect on drought. The SPEI 
uses the "climatic water balance" component, which is the difference between precipitation and 
potential evapotranspiration rather than precipitation as the input variable. Specifically, the climatic 
water balance compares water availability with atmospheric evaporative demand and therefore 
provides reliable drought severity measurements rather than only considering precipitation as the 
input variable (Begueria et al., 2014). The output of the SPEI is applicable for all climate regimes, 
and the results are comparable because they are standardized (WMO and GWP, 2016). The SPEI 
can be calculated at different time scales. Therefore, it is possible to detect both short- and 
long-term drought (Vicente-Serrano et al., 2010). The calculation of the SPEI is based on monthly 
precipitation and monthly potential evapotranspiration. 

The SPEI package version 1.7 in R software was used to calculate the monthly SPEI values in 
the study area. The negative values of the SPEI indicate drought conditions. The lower the negative 
values, the worse the drought conditions. The multiscale features of time for the SPEI values at a 
3-month time scale (May for spring, August for summer, November for autumn, and February for 
winter in our case) provide drought in the seasonal context. Similarly, the SPEI values at the 
12-month time scale (December SPEI of each year) of each year indicate the whole-year water 
conditions caused by drought. The 3-month and 12-month time scales reflect the agricultural and 
hydrological drought conditions, respectively. The 3-month time scale also describes the short- and 
medium-term soil moisture and the seasonal precipitation estimation. Similarly, the 12-month time 
scale observes the long-term precipitation pattern over a period. The 12 consecutive months can 
also provide annual moisture conditions (WMO, 2012; Wang et al., 2014). This study examined the 
annual and seasonal drought characteristics in the KRB and focused on moderate, severe, and 
extreme drought events based on the SPEI values illustrated in Table 3 (McKee et al., 1993; 
Vicente-Serrano et al., 2010). 


Table 3 Classification of drought based on the standardized precipitation evapotranspiration index (SPEI) values 


SPEI value Category 
> -0.99 Near normal 
—1.49- —1.00 Moderate drought 
—1.99- —1.50 Severe drought 
< -2.00 Extreme drought 


Based on this classification, we extracted 3-month (seasonal) and 12-month (annual) time scales 
of the monthly SPEI values from 1987 to 2017. The livelihood in the KRB is based on agriculture 
(Neupane et al., 2013); thus, studying drought at the seasonal scale can be essential (Tan et al., 
2015). Therefore, the monthly SPEI values at the 3-month time scale were used to calculate the 
drought event characteristics (intensity, duration, severity, and frequency). The annual SPEI was 
indicated by the December SPEI of each year. Similarly, the seasonal SPEI was obtained in May 
for spring, August for summer, November for autumn, and February for winter. 


2.3.2 Estimation of potential evapotranspiration 


To calculate the SPEI values, it is necessary to estimate potential evapotranspiration from the 
available climate data. There are several methods to estimate this variable, and some commonly 
used methods are the Thornthwaite method (Thornthwaite, 1948), Hargreaves method (Hargreaves 
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and Samani, 1985), and Penman-Monteith method (Allen et al., 1994; Allen et al., 1998). Each 
method has different data inputs. Several studies have used these methods based on input data 
availability; for example, Tan et al. (2015) and Mohsenipour et al. (2018) used the Thornthwaite 
method, Prabnakom et al. (2018) and Bisht et al. (2019) used the Hargreaves method, whereas Gao 
et al. (2017) and Huang et al. (2018) used the Penman-Monteith method. Among these methods, 
the Penman-Monteith method (Penman-Monteith FAO-56) is recommended by the Food and 
Agriculture Organization (FAO) because it has a greater potential for estimating accurate PET in a 
wide range of climate conditions (Allen et al., 1998). However, this method requires a number of 
inputs, such as wind speed, net radiation, temperature, and relative humidity, which are not readily 
available, especially in developing countries across the globe (Farmer et al., 2011). Therefore, the 
alternative methods that require relatively few data are the Thornthwaite and Hargreaves methods; 
they are widely used for the estimation of potential evapotranspiration (Xu and Singh, 2001; Tam et 
al., 2019). 

Thomthwaite and Hargreaves methods are the temperature-based approaches for calculating 
potential evapotranspiration. The Thornthwaite method requires the average temperature and 
latitude (calculated sunshine duration) to calculate potential evapotranspiration. The Hargreaves 
method uses the minimum and maximum temperatures and extraterrestrial radiation as inputs 
(Hargreaves and Samani, 1985). This method is recommended as the standard method for 
computing potential evapotranspiration when only temperature data are available (Bandyopadhyay 
et al., 2009; Begueria et al., 2014; Khan et al., 2020). It is also recommended by the FAO as an 
alternative method (Allen et al., 1998), which is comparable to the findings of the 
Penman-Monteith method (Droogers and Allen, 2002). Stagge et al. (2014) studied the sensitivity 
of the SPEI to potential evapotranspiration estimation methods. Their study includes five potential 
evapotranspiration estimation methods, namely, the Thornthwaite, Hargreaves, Penman-Monteith 
Hargreaves, Priestley Taylor, and Penman-Monteith FAO-56. Among these, the SPEI estimated 
from the Hargreaves and the other three methods (Penman-Monteith Hargreaves, Priestley Taylor, 
and Penman-Monteith FAO-56) are more internally consistent than that estimated from the 
Thornthwaite method. Therefore, the Hargreaves method is the most useful and consistent method 
to estimate potential evapotranspiration in data-scarce areas (Stagge et al., 2014). 

Therefore, our study also adopted the Hargreaves method to estimate potential 
evapotranspiration because of data limitations. The Hargreaves method has also been successfully 
applied in Nepal. For example, Bhatt et al. (2014) estimated potential evapotranspiration of the 
KRB; Dahal et al. (2016) studied the impact of climate change on water availability in the Bagmati 
River Basin; Penton et al. (2016) used this method for modeling the hydroclimate and stream flow 
in the KRB; Aadhar and Mishra (2017) studied real-time drought analysis and monitoring in South 
Asia (including Nepal); and Hamal et al. (2020) assessed drought impacts on crop yields in Nepal. 
The equation of the Hargreaves method can be expressed as follows: 


ET, = 0.0023 x RA(T. + 17.8) xTD (1) 


avg 


where ETo is the mean monthly reference or potential evapotranspiration (mm/month); RA is the 
extraterrestrial radiation (MJ/(m?-month)) estimated from the latitude and month of the year; Tavg is 
the average monthly temperature (°C); TD is the mean monthly temperature range (°C) (difference 
between mean maximum and mean minimum temperatures); and 0.0023 is the coefficient proposed 
by Hargreaves and Samani (1985). 
2.3.3 Definition of drought characteristics 
A drought event is defined as a period in which the drought index value is continuously negative 
and reaches the value of —1.00 or less. Drought starts when the value falls below zero and ends 
when it reaches a positive value following the value of —1.00 or less (McKee et al., 1993). 

Drought frequency is defined as the number of drought events to the total years of the study 
period. It is calculated as follows: 


F, =>. x100%, (2) 
N 


S 
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where F, is the drought frequency (%); ns is the number of drought events; N; is the total years for 
the study period; and s is the station. 

Drought duration is the number of months between the start and end of a drought event. The start 
month will be included and the end month will be excluded during the counting (Spinoni et al., 
2014). Severity is the sum of all SPEI values (absolute) during a drought event; it is defined as 
follows: 


S, = Y` SPEI, 


j=l 


(3) 


e 


where e is the drought event; j is the month; and m and Se are the duration and severity of a drought 
event e, respectively. 

Drought intensity (DIe) is obtained by dividing severity by duration (Eq. 4). The larger the 
intensity, the greater the drought severity. 


DI, ==. (4) 
m 
The drought station proportion is the number of drought stations on a time scale relative to the 
total number of stations studied. It provides the spatial extent of drought events during a period 
over the stations (Gumus and Algin, 2017). It is calculated as follows: 


n. 
P, = x100%, (5) 
N; 


where Pj; is the drought station proportion at a given time scale (%); j is the time scale; njis the 
number of droughts in the given time scale; and Nj is the total number of stations. 

2.3.4 Trend analysis 

The M-K test and Sen's slope methods (Salmi et al., 2002) were used for trend analysis. These 
methods are widely used in detecting trends in the hydrometeorological data series (Burn and 
Elnur, 2002; Tan et al., 2015), including drought index series (Damberg and AghaKouchak, 2014). 
The M-K test is a nonparametric test and does not require any specific data distribution. The Sen's 
slope method is also not highly affected by single data errors or outliers (Salmi et al., 2002). These 
tools can help to obtain the magnitude and trend significance in the time series data. The 
significance of the monotonic positive or negative trend was determined by the M-K test, and the 
Sen's slope provided the magnitude of change. The positive and negative signs in the value indicate 
increasing and decreasing magnitudes, respectively. The M-K test uses the hypothesis testing 
approach towards the monotonic trend in the time series data. It tests the null hypothesis (Ho) if a 
trend is not observed in the time series data as well as the alternative hypothesis if a trend is not 
observed in the data series at the 0.001, 0.01, and 0.05 levels of significance. The M-K test statistic 
S is calculated as follows: 


S = 2 sgn(x, -x,) (6) 


where x; and x; are the annual values in years j and k (j>k), respectively. 
1, if x,- x, >0 
sgn(x, — X,) =40, if x, -—x, =0. (7) 
-1, if x,-x<0 
When the number of values in data series (n) is less than or equal to 9, the absolute S value is 


compared directly to the theoretical distribution of S derived by the M-K method (Gilbert, 1987). 
When n=>10, the normal approximation test is used, and the variance of S is calculated as follows: 
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VAR(S) == n(n—1)(2n+5) Yt, (t, -1)(2e, +5) , (8) 


p=1 


where VAR(S) is the variance of S; q is the number of connected groups in series; and tpis the 
number of data values in the pè group. The test statistic Z is calculated using S and VAR(S) as 
shown below: 


S-1 
VARS)’ 
Z=40, ifS=0 . (9) 
S+1 


JVAR(S) 


The Z value is used to evaluate the presence of a statistically significant trend. Positive and 
negative Z values indicate the increasing and decreasing trends, respectively. The increasing and 
decreasing trends (a two-tailed test) are tested at the a level of significance. Hois rejected if the Z 
value is greater than the Zi-«2 value, where Zı-a2is obtained from a standard normal cumulative 
distribution. In this study, we tested the significance at 0.001, 0.01, and 0.05 levels. 

The trend magnitude (Q) is quantified by the Sen's slope method and computed as follows: 


ifS>0 


if S<0 


X,—X 
Q, E (k+ j; i=1, 2,3, N), (10) 
j- 
where j and k are the data series that have n number of values, and x; and xxare the data values of 
these series. The slope observations (N) are calculated as N=n(n—1)/2. The median of N values of 
slope observations gives Sen's estimator of slope Q;. The N values of Q; are ranked from ascending 
to descending order, and then the Sen's estimator is expressed as follows: 


se 1 oe 
Q= Quay if Nisodd or Q =5 (Qw + Qsan) ) if Niseven . (11) 


Finally, Qis tested by a two-sided hypothesis test at a 100% (1-a) confidence interval, and the 
true slope is obtained. If the slope is significant at the 95% confidence interval by means of a 
nonparametric technique based on a normal distribution, then the quantity of a (C,) is computed as 


follows: 
C, = Zian y VAR(S) , (12) 


where VAR(S) is defined in Equation 8 and Zı-wa2 is estimated from the standard normal 
distribution. 
Z represents the standard normal deviate and a is set as 0.05, then the lower and upper 
confidence limits (Mı and M2, respectively) are calculated as follows: 
M,=(N-C,)/2, (13) 
M, =(N+C,)/2. (14) 
The lower and upper limits of the confidence interval, Qmin and Qmax, are the Mi" largest and the 
(M2+1)" largest of the N ordered slope estimates Qi, respectively. If Mı and M are not whole 
numbers, then the lower and upper limits are interpolated, respectively (Salmi et al., 2002). In this 
study, we used the Microsoft Excel program MAKESENS 1.0 developed by the Finnish 
Meteorological Institute in 2002 for the analysis (Salmi et al., 2002). 


3 Results 


3.1 Spatial and temporal variations in the SPEI 


The annual variation in the SPEI showed a decreasing trend at the regional and basin levels of the 
KRB from 1987 to 2017 (Fig. 2). The negative SPEI values represent the drought conditions (as 
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shown in Table 3). The temporal variation from the SPEI analysis results revealed that the 
mountain region experienced moderate drought in 1992, 2011, 2014, and 2015, whereas severe 
drought occurrences were observed in 2005 and 2009. Similarly, in the hilly region, moderate 
drought was observed in 1992, 2008, 2009, 2015, and 2017, and severe drought occurred in 2012. 
Moderate drought was observed in 2009 and 2015 and severe drought was observed in 1992, 1994, 
and 2012 in the Terai region. 


—e— SPEI ---- Moderate drought ---- Severe drought ---- Extreme drought 
35 3 
(a) Mountains Trend= —0.30/10a (b) Hills Trend= —0.60/10a 
24 P>0.05 2 ; P<0.05 


SPE] value 


3 T T T T T T T 1 3 T T T T T T 1 
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24 P>0.05 2 P>0.05 
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Fig. 2 Annual variation of the standardized precipitation evapotranspiration index (SPEI) values at the 12-month 
time scale at the regional (a—c) and basin (d) levels of the KRB from 1987 to 2017. The trend of the SPEI per 
decade (Mann-Kendall test) and the statistical level of significance are also shown. 


Spatially, annual variation in moderate drought was observed in the mountains, hills and Terai 
during the recent year (e.g., 2015). In 1992, drought was moderate in the mountains and hills, 
whereas it was severe in the Terai. In 2009, drought was severe in the mountains and moderate in 
the hills and Terai. Similarly, the analysis showed severe drought in the hills and Terai in 2012. Our 
results revealed that the drought conditions in the mountain, hill and Terai regions are independent. 
However, for certain years, the occurrence of drought might be related, which was not consistent 
over the study period of 1987-2017. 

The overall results indicated that drought occurrence was frequent after 2000. Moderate drought 
occurrence was relatively higher in the hill region, and severe drought was higher in the Terai 
region. Spatially, the consistency in annual drought occurrence was higher within the hills and Terai 
among the three regions. The annual SPEI values showed a significant decreasing trend every ten 
years in the hills, indicating the great possibility of increasing drought tendency in the future. 

Similar to the annual trend, the seasonal variation in the SPEI also showed a decreasing trend in 
the regional and basin levels of the KRB from 1987 to 2017 (Fig. 3). In the mountain region, 
extreme drought was observed in autumn (2009); severe drought was observed in spring (1992 and 
2011), summer (1993 and 2009), and winter (1999 and 2006); and moderate drought was also 
noticed in spring (2005, 2007, 2014, and 2016), summer (1992, 2005, and 2011), autumn (1994, 
2000, 2002, 2005, and 2015) and winter (2008, 2009, 2014, and 2017). 

In the hilly region, severe drought occurred in spring (1996), summer (2006 and 2017), autumn 
(2000), and winter (2006 and 2009); moderate drought also occurred in spring (1992, 1995, 1999, 
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2012, and 2014), summer (2008, 2009, and 2012), autumn (1988, 1997, 2008, 2009, 2012, 2015, and 
2017) and winter (1999, 2001, and 2017). 


---- Moderate drought ---- Severe drought ---- Extreme drought 
3 C3 Spring SPEI WE Summer SPEI CJ Autumn SPEI WEE Winter SPEI 
> (a) Mountains . 347 (b) Mountains 3 
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Fig. 3 Seasonal variation in the SPEI values at the 3-month time scale at the regional (a-f) and basin (g and h) 
levels of the KRB from 1987 to 2017. The trend of the SPEI per decade (Mann-Kendall test) and the statistical level 
of significance are also shown. NS, not significant. 


In the Terai region, an extreme autumn drought was observed in 2009. Similarly, severe drought 
was observed in spring (1995 and 1996), summer (1992, 1994, and 2012), autumn (2005), and 
winter (2006, 2017); and moderate drought was also noticed in spring (1991, 1992, 1994, and 2012), 
summer (2006, 2013, and 2015), autumn (2008, 2012, and 2015), and winter (1999 and 2001). The 
SPEI values showed a decreasing trend in all four seasons during the entire study period. However, 
the trend is only significant during the summer season in the hill region and during the winter season 
in the mountain and hill regions. 

Spatially, the seasonal drought showed variations in the study period. In some years, a similar 
seasonal drought pattern was observed for all the regions. For instance, in the mountain, hill and 
Terai regions, severe winter drought was observed in 2006 and moderate autumn drought was 
observed in 2012. Similarly, in some years, although all regions showed seasonal drought, the 
pattern was consistent in at least two regions. For example, in 1992, moderate spring drought was 
observed in the hills and Terai while severe drought was observed in the mountains. A similar pattern 
was also observed for winter drought in 1999. Further, in 2009, an extreme autumn drought was 
observed in the mountains and Terai while moderate drought was observed in the hills. In 2017, 
moderate winter drought was observed in the mountain and hill regions, while severe drought was 
observed in the Terai region. 

Furthermore, the results showed seasonal drought in at least two regions. The drought severity 
was similar in some cases (mostly in the hill and Terai regions) but different in others. For example, 
the hill and Terai regions showed severe spring drought in 1996, moderate winter drought in 2001, 
moderate autumn drought in 2008 and 2015, and moderate spring drought in 2012. In the same way, 


JOURNAL OF ARID LAND 


moderate summer drought was observed in the mountain and hill regions in 2011 and moderate 
spring drought was observed in 2014. 

In some years, different severities of seasonal drought were also observed in two regions. For 
instance, during the autumn of 2000, moderate drought was observed in the mountains while severe 
drought was observed in the hills. During the same season (autumn) in 2005, the mountains showed 
moderate drought and the Terai showed severe drought. Similarly, in 1995, moderate spring drought 
was detected in the mountain region and severe drought was detected in the Terai region. During the 
summer of 1992, moderate drought and severe drought were observed in the mountain and Terai 
regions, respectively. In addition, during the summers of 2006 and 2009, severe drought was 
observed in the hills and mountains while moderate drought was observed in the Terai and hills. 
Similarly, during the winter of 2009, moderate drought was observed in the mountains and severe 
drought was observed in the hills. 

The overall results of the seasonal variation in the SPEI indicated that both moderate drought 
(most common) and severe drought occurred in all regions of the KRB. In some years, a higher 
spatial consistency in seasonal drought was observed in two regions than in all regions. A significant 
decreasing trend in the SPEI was observed in the hill region in summer but in the mountain and hill 
regions in winter. In the entire KRB region, a significant decreasing trend of seasonal SPEI values 
was observed in the summer and winter seasons. 


3.2 Analysis of drought characteristics 


3.2.1 Duration and severity of drought events 


The average drought duration calculated from the SPEI for the KRB during the period 1987—2017 
was 5.8 months. In this basin, the maximum drought duration (for approximately 15.0 months) was 
observed during the period 2009-2010. Spatially, substantial differences in the average drought 
duration from 1987 to 2017 were not observed between the mountain and hill regions, although the 
values were relatively lower (by a month) in the Terai region. The maximum drought duration 
(from 2007 to 2009) was greatest in the hill region. The occurrence of the maximum drought 
duration in the mountain region was observed twice from 2005 to 2013, whereas it was only 
observed in the mid-1990s in the Terai region (Table 4). 

The average and maximum severity values of drought events in the KRB from 1987 to 2017 
were 5.80 and 16.05, respectively. Spatially, the average severity of drought events was highest in 
the hilly region and lowest in the Terai region, although the magnitude of the maximum severity 
was highest in the Terai region and lowest in the mountain region. The longest duration of the 
maximum severity was observed in the mid-1990s in the Terai region but in the mid-2000s and late 
2000s in the mountains and hills. 


Table 4 Summary of the duration, severity, and intensity of drought events identified from the SPEI values at the 
3-month time scale from 1987 to 2017 at the regional and basin levels of the KRB 


Maximum drought 


Region pnie Starting—ending time Severity Starting-ending time Intensity Starting-ending time 
Mountain 14.0 pA TEEN Sah 15.49 Apr 2005-May 2006 1.74 Jun 2009-Nov 2009 
Hill 17.0 Dec 2007—Apr 2009 15.72 Dec 2007—Apr 2009 1.39 Jan 2006—Mar 2006 
Terai 16.0 Apr 1994—Jul 1995 17.28 Apr 1994—Jul 1995 1.84 Nov 2005 
KRB 15.0 Jan 2009—Aug 2010 16.05 Feb 2012—Mar 2013 1.59 Jan 2006—Mar 2006 
F Average drought Drought frequency 
Region Duration Severity Intensity (%) 
(month) 
Mountain 5.7 5.37 0.94 74 
Hill 5.9 5.60 0.94 77 
Terai 4.7 4.77 1.01 81 


KRB 5.8 5.80 1.00 71 
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3.2.2 Intensity and frequency of drought events 

The average and maximum drought intensities calculated from the SPEI in the KRB during the 
period 1987-2017 were 1.00 and 1.59, respectively. The maximum intensity was observed over a 
3-month duration in 2006. Spatially, the Terai region showed the highest magnitude of the average 
and maximum drought intensity compared with the mountain and hill regions. The highest drought 
intensity in the Terai region was observed in a single month duration (November) of 2005. The 
magnitude of the average drought intensity in the mountain and hill regions was 0.94, while the 
maximum intensity was higher in the mountain region (1.74). The duration of the maximum 
drought intensity in the mountain, hill, and Terai regions was six, three, and one months, 
respectively. The frequency of drought events in the KRB from 1987 to 2017 was 71%. Spatially, 
the frequency of drought events was 74% in the mountain region, 77% in the hill region, and 81% 
in the Terai region, which indicated that there was a general increment of drought frequency events 
from the mountain region to the Terai region. 

The overall results from the SPEI analysis indicated that the maximum duration, intensity, and 
severity of drought events were mainly observed after 2000 in the KRB region. Spatially, the onset 
of the maximum duration and severity of drought events occurred either in the winter or spring 
season in all regions, including the entire KRB, and ended mostly in the spring season of the next 
year. Similarly, the maximum drought intensity was observed in summer and late autumn in the 
mountain region, in the mid-winter and early spring in the hill region and the entire KRB, and in 
the late autumn in the Terai region. 


3.3 Variation in drought station proportion 


The monthly drought station proportion was calculated from the monthly SPEI values using 
Equation 5 for the period 1987—2017. The annual and seasonal drought station proportions for the 
same period were calculated by averaging the monthly proportions of a year and their respective 
seasons. The M-K trends of these proportions at the monthly, interannual, and interseasonal scales 
from 1987 to 2017 are illustrated in Table 5. In the mountain region, the proportion at the monthly 
scale did not show trends in any month. However, at the interseasonal and interannual scales, 
significant increasing trends of 4.1%/10a and 5.8%/10a (P<0.05) were found, respectively. 
Similarly, in the hill region, a significant increasing trend was found in January (5.3%/10a; 
P<0.05), March (3.0%/10a; P<0.05), July (5.1%/10a; P<0.05), August (7.6%/10a; P<0.05), and 
September (7.0%/10a; P<0.01); no trend was found in other months. At the interseasonal scale, the 
proportion showed an increasing trend in all seasons but was only significant in summer 
(7.2%/10a; P<0.05), autumn (4.5%/10a; P<0.05), and winter (5.7%/10a; P<0.05). At the 
interannual scale, the proportion also showed a significant increasing trend (6.7%/10a; P<0.01). 
Similarly, in the Terai region, the proportion showed a significant decreasing trend (—4.8%/10a; 
P<0.05) in June and no trends were observed in other months. At the interseasonal scale, a positive 
trend was observed in spring and autumn, a negative trend was observed in winter, and no trend 
was observed in summer. Similarly, at the interannual scale, no trend was observed. 

The overall results indicated that the spatial extent of drought has a significant increasing trend 
in summer, autumn, and winter in the hill region. In the entire KRB region, the drought station 
proportion calculated from the SPEI showed an increasing trend in January and February and from 
July to November; no trends were observed from March to June and in December. However, only 
the trends in February, August, and September were statistically significant (P<0.05; Table 5). At 
the interseasonal scale, the drought station proportion in all seasons showed an increasing trend but 
was not statistically significant. However, the annual drought station proportion showed a 
significant increasing trend (4.2%/10a; P<0.05) from 1987 to 2017. 


4 Discussion 


Our study showed that the drought and its characteristics exhibited significant spatial and temporal 
variations in the KRB from 1987 to 2017. These variations are directly and indirectly influenced by 
many factors, such as precipitation, temperature, and topography. The variations on these factors, 
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especially the decreasing precipitation and increasing temperature, are the contributing factors of 
drought events. Thus, the variations also lead to differences in the drought intensity and duration 
(Trenberth, 2011; Wang et al., 2019). A global-level study on drought revealed that since the 
mid-1980s, the drought tendency has been increasing because of the increasing temperature trend 
through evapotranspiration mechanism (Dai, 2011; Vicente-Serrano et al., 2011). Similarly, the 
reduced precipitation frequency triggers and enhances the drought (Potop et al., 2014). The 
topographic diversity in a region can also lead to drought. Topographic features can create a rain 
shadow region that receives low precipitation and reinforce the drought (Dixit et al., 2009). Our 
results exhibited more significant increasing trend of annual and seasonal temperatures in the hill 
region than in the other regions (Table 6). The winter precipitation showed significant decreasing 
trends in all regions while the annual trend was only observed in the hills. To better understand the 
drought conditions at the spatial and temporal scales, we discussed some influencing factors and 
their implications as follows. 


Table 5 Mann-Kendall (M-K) trends in drought station proportion series calculated from the SPEI values at the 
3-month time scale at the regional and basin levels of the KRB from 1987 to 2017 


M-K trend in drought station proportion (%/10a) 


Time 
Mountains Hills Terai KRB 
January 0.0 5.3" 0.0 2.2 
February 0.0 0.0 0.0 2.9" 
March 0.0 3.0" 0.0 0.0 
April 0.0 0.0 0.0 0.0 
May 0.0 0.0 0.0 0.0 
June 0.0 0.0 4,8" 0.0 
July 0.0 5.1° 0.0 3.3 
August 0.0 7.6" 0.0 5.7° 
September 0.0 7.0" 0.0 4.8" 
October 0.0 0.0 0.0 3:3: 
November 0.0 0.0 0.0 2.9 
December 0.0 0.0 0.0 0.0 
Spring 0.0 3.0 1.6 0.7 
Summer 4.1 7.2" 0.0 3.8 
Autumn 0.0 4.5" -2.5 4.6 
Winter 0.0 5.7° -1.4 3.9 
Annual 5.8° 6.7" 0.0 4.2" 


Note: ™ and * represent M-K trends statistically significant at 0.01 and 0.05 levels, respectively. 


Table 6 M-K trends of seasonal and annual precipitation and temperature at the regional and basin levels of the 
KRB from 1987 to 2017 


Season 
Region Climate variable Annual 
Spring Summer Autumn Winter 
Precipitation (mm/a) —6.40 —7.70 4.80 -11.70° —25.10 
Mountain 

Temperature (°C/a) 0.20 0.10 0.10 0.30 0.20 

Precipitation (mm/a) 1.70 68.30 —0.40 -10.90° 84.40" 
Hill 

Temperature (°C/a) 0.50™ 0.40™ 0.30™ 0.50% 0.40" 

Precipitation (mm/a) 11.80 —59.60 -13.40 —8.00° —59.00 
Terai 

Temperature (°C/a) 0.20 0.20™ 0.10 0.03 0.10 

Precipitation (mm/a) 2.20 —64.00 —12.50 —10.50° —60.60 
KRB 

Temperature (°C/a) 0.30™ 0.30°" 0.20% 0.20* 0.30% 

Note: “**, ™*, and * denote significance levels of 0.001, 0.01, and 0.05, respectively. 
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4.1 Spatial and temporal drought trend variations in the KRB 


One of the purposes of this study is to scrutinize the spatial and temporal variations in drought 
based on the SPEI. Our findings revealed that the drought trend in the basin is increasing spatially 
and temporally. Significant increasing drought trends were observed at the annual and seasonal 
(summer and winter) scales in the hills and in winter in the mountains. This increasing drought 
trend indicates that the occurrence of drought is likely to be more frequent and thus can negatively 
impact water resources, agriculture, and livelihood in the KRB. The drought findings in the hill 
region resembles the findings of some quantitative and qualitative studies conducted in the hill 
region of Nepal. For instance, Ghimire et al. (2010) identified that 63% of farms in the hills of 
Nepal are highly vulnerable to drought; Kafle (2014) identified incidents of extreme drought events 
from five meteorological stations in the mid- and far-western hill regions of Nepal; Sapkota et al. 
(2016) mentioned the severe drought-affected area in Kavre District in the hilly region of the KRB 
and explored climate change vulnerability; Adhikari (2018) studied the occurrence of drought and 
its impact on rainfed agriculture in the western hills of Nepal. 

Several studies have explained that increasing temperature and decreasing precipitation are 
directly linked to drought variations (Challinor et al., 2009; Sigdel and Ikeda, 2010; Liu et al., 
2015; Dehghan, 2020). The variation in these climatic variables creates higher atmospheric 
evaporative demand and long dry spell durations, thereby increasing the dynamics and magnitude 
of drought (NCVST, 2009; Vicente-Serrano et al., 2014; Wang et al., 2014). Based on these facts, 
the annual and seasonal precipitation and temperature trends of the three topographic regions 
(mountain, hill, and Terai) and the entire KRB from 1987 to 2017 are shown in Table 6. The 
analysis in the entire KRB in our study is consistent with the findings of Baniya et al. (2018) that 
the annual precipitation decreased at the rate of 5.12 mm/a and the temperature increased at the rate 
of 0.03°C/a from 1982 to 2015 in Nepal. The increasing and decreasing trends in temperature and 
precipitation at both annual and seasonal scale also fit with the respective SPEI trends from our 
findings (Figs. 2 and 3), which indicate a high drought risk. The other causes of variations in 
drought might be the geographical features of the KRB. As the elevation increases from the Terai 
(south) to the mountains (north), the climate varies from tropical to warm temperate (Shrestha et 
al., 2017). The South Asian monsoon influences the climate of the southern and northern parts, 
which lie in rain shadow areas. The precipitation pattern in the basin is linked with summer 
monsoon-southeasterly and winter monsoon-westerly circulations (Sigdel and Ikeda, 2010). 

Due to substantial differences in topography, the spatial and temporal variations in precipitation 
are observed within a short distance. In the mountain region, even after the monsoon, the 
north-facing slopes have a relatively higher humidity than the south-facing slopes because of 
diffuse solar radiation in the north-facing slopes. On a hill, a sudden cloudburst can cause intense 
precipitation within a day, whereas in the rain shadow areas of the Tibetan Plateau, the conditions 
may be dry (Dixit et al., 2009). Therefore, local geographical features, orography, and complex 
terrain produce varied microclimates that might cause drought conditions in the basin (Kansakar et 
al., 2004; Sigdel and Ikeda, 2010; Zuo et al., 2018). 


4.2 Variation of drought characteristics in the KRB 


Another key purpose of this study was to detect the variation of drought characteristics in the KRB 
using the SPEI values. The analysis of the SPEI values showed the highest frequency and intensity 
(average and maximum) of drought in the Terai region. Similarly, the highest duration of the 
maximum drought was observed in the hill region. Saravi et al. (2009) indicated that a drought 
event with the longest duration may not necessarily be intense enough to result in the greatest 
impacts. In our case, the Terai region had only a 1-month drought duration in 2005 with the 
maximum drought intensity of 1.84, whereas the mountain and hill regions had a relatively longer 
drought duration (Table 4). 

Our study found that the magnitude of drought intensity was highest in the Terai region, which 
might have impacted the growth of winter crops, especially wheat. Our finding is similar to the 
results of Wu et al. (2019) that severe drought occurred in the downstream of the KRB from 2000 
to 2014. The study by Wu et al. (2019) in the transboundary regions of the KRB revealed that the 


JOURNAL OF ARID LAND 


highest drought intensity is directly related to the higher drought severity and thus it has a 
significant impact on crops and water facilities. Our study observed the occurrence of drought in 
different seasons, especially winter and spring in the hills and summer, winter, and spring in the 
mountains, Terai, and the entire KRB. These seasons are highly important for producing major 
cereals crops, such as wheat, maize, and rice. Therefore, any kind of drought conditions might 
affect production. These findings support the results of Joshi (2018) mentioned by Adhikari (2018), 
who revealed that the eastern and central regions of Nepal faced droughts in 1992, 1994, 1997, 
2002, 2008, 2009, 2012, and 2013 because of inadequate and poor distribution of precipitation, 
which caused major crop loss of approximately 3.85x10° t. The seasonal occurrence of drought is 
also consistent with the findings of Baniya et al. (2019) that during the period of 1982-2015, 
several severe and moderate droughts occurred during the decades of the 1990s, 2000s, and 2010s 
at the annual and seasonal scales. Similarly, a study in Nepal also showed moderate to extreme 
drought intensity from 1972 to 2003 (Sigdel and Ikeda, 2010). Some other findings also revealed 
that the western part of Nepal is also experiencing extreme drought conditions, especially in 
summer and winter (Wang et al., 2013; Kafle, 2014). A study using the SPI values from 1981 to 
2012 in the Gandaki River Basin of Nepal showed that 2004, 2005, 2006, and 2009 were the years 
with severe drought (Dahal et al., 2016). These results on drought from several studies indicated 
that different regions and river basins of Nepal along with the KRB experienced drought from the 
1980s to 2010s. 


4.3 Implications of drought variations in the KRB 


The implications of drought occurrence and variations in any location at a given time can be 
determined by the climate features and biophysical conditions of the location (Wilhite, 2000). 
Therefore, these conditions can also evaluate the seasonal and annual drought in the KRB. The 
3-month time scale highlights the soil moisture conditions and the 12-month time scale reflects the 
water levels in the source regions; these time scales are most effective in the agricultural region. 
Our findings revealed that the annual and seasonal drought trends were increased in the KRB from 
1987 to 2017 (Figs. 2 and 3), which can reduce the soil moisture and decrease the water source 
levels and water availability (Wilhite, 2000; Neupane et al., 2013; Amrit et al., 2018). It is accepted 
that the increased intensity, duration, and frequency of drought can affect water sources and 
agricultural productivity (Peterson et al., 2012; Duan and Mei, 2014). A review study by Adhikari 
(2018) on the hills of Nepal showed that the climate-induced drought profoundly affects 
agricultural systems. The KRB region has a significant quantity of arable land and therefore is 
more sensitive to drought. Aryal (2012) mentioned that crops are sensitive to the water 
requirements during their lifecycle. In the stage of crop growth and development, crops have a 
high-water demand and any delay of water supply would affect the growth of crops and decrease 
the production potential (Paudyal et al., 2001). A study conducted by Bhatt et al. (2014) mentioned 
that rice, maize, and wheat cultivated below 1100, 1350, and 1700 m a.s.l. in the KRB, 
respectively, suffer from drought caused by increasing temperatures during the flowering and yield 
formation phases. Thus, the delay in fulfilling water demand at the growth stage will ultimately 
affect the yield. Malla (2008) revealed that in Nepal, the number of rainy days is decreasing and the 
rainy days are highly intense followed by prolonged drought. In the KRB region, the sectoral water 
demand increases from the mountain region to the Terai region; it is highest in the agricultural 
sector (92%), followed by the domestic (6%) and industrial (2%) sectors (Chinnasamy et al., 2015). 
Temporally, the water demand is higher in the dry season than in the wet season. In this scenario, 
our results imply that drought may decrease the water availability for all sectors in the future; 
therefore, necessary policies and plans for water management, such as increasing water storage 
capacity and other drought mitigation strategies, should be adopted. 


5 Conclusions 


In this study, the spatial and temporal characteristics of drought in the KRB were assessed from 
1987 to 2017. The spatial and temporal variations in the SPEI values indicated that drought 
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occurrence was more frequent after 2000. The moderate and severe drought occurrences were 
relatively higher in the hill and Terai regions, respectively. The annual SPEI values showed a 
significant decreasing trend per decade in the hill region and the entire KRB, implying an 
increasing drought tendency in the future, which could be further exacerbated due to climate 
change. The rank of the magnitude of the drought trend was hills>mountains>Terai. The seasonal 
variation in the SPEI indicated that both moderate (most common) and severe droughts occurred in 
all regions of the KRB. These droughts had a significant increasing trend in winter over the 
mountain and hill regions and in summer and winter over the hill region. The overall rank of the 
magnitude of the spring drought trend was Terai>mountains>hills; rank of summer drought trend 
was hills>mountains>Terai; rank of winter drought trend was mountains>hills>Terai; and rank of 
autumn drought trend was the hills and Terai>mountains. Therefore, the regions that are more 
likely to be affected by seasonal drought should focus on essential mitigation measures. 

The analysis of the drought characteristics indicates that the maximum duration, intensity, and 
severity of drought events were observed in the KRB after 2000. These drought characteristics 
mainly appeared either in winter or spring of the current year and mostly ended in spring of the 
next year in the entire region of the KRB. The Terai region had the highest drought frequency 
(81%) and average drought intensity (1.01), while the hill region had the highest average drought 
severity (5.60) and drought duration (5.9 months). Similarly, the hill region had the longest 
maximum drought duration (17.0 months) and the Terai region had the maximum drought severity 
(17.28). The drought station proportion calculated from the SPEI values at a 3-month time scale 
showed the highest significant increasing trend in the hill region at the monthly (7.6%/10a in 
August), seasonal (7.2%/10a in summer), and annual (6.7%/10a) scales. The spatial extent of 
drought decreased only in the Terai in June in summer (P<0.05), autumn and winter during the 
entire study period. 

In this study, we attempted to observe drought conditions in different ecological regions of the 
KRB using the SPEI values. The overall findings revealed that drought is significant in the KRB; 
the results provide a basis for long-term monitoring of drought in the future. This study could 
provide resource planners and policy-makers information on drought conditions under a changing 
climate, which can help them make informed climate decisions and implement them. Moreover, a 
further detailed study should be conducted to understand the drought impact at the district level. In 
addition, further research is necessary to investigate the detailed impacts of increasing drought 
trends during crop growth phases and water demand in scenarios of population growth. 
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